Non-magnetic and magnetic thiolate-protected Au25 superatoms on Cu(lll), Ag(lll) 

and Au(lll) surfaces 
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Geometry, electronic structure, and magnetic properties of methylthiolate-stabilized AU25L18 and 
MnAu24Li8 (L = SCH3) clusters adsorbed on noble-metal (111) surfaces have been investigated by 
using spin-polarized density functional theory computations. The interaction between the cluster 
and the surface is found to be mediated by charge transfer mainly from or into the ligand mono- 
layer. The electronic properties of the 13-atom metal core remain in all cases rather undisturbed as 
compared to the isolated clusters in gas phase. The Au2sLig cluster retains a clear HOMO - LUMO 
energy gap in the range of 0.7 eV to 1.0 eV depending on the surface. The ligand layer is able to 
decouple the electronic structure of the magnetic MnAu24Lig cluster from Au(lll) surface, retaning 
a high local spin moment of close to 5 fiB arising from the spin-polarized Mn(3d) electrons. These 
computations imply that the thiolate monolayer-protected gold clusters may be used as promising 
building blocks with tunable energy gaps, tunneling barriers, and magnetic moments for applications 
in the area of electron and/or spin transport. 

PACS numbers: 73.22.-f, 36.40.-c, 61.46.-w 



Ligand monolayer protected gold clusters (MPCs) have 
received a lot of interest since mid-1990's due in part 
to their enhanced stability, easy synthesis, rather well- 
defined monodispersity and potential applications in nan- 
otcchnology. For example, the MPCs may be used as cat- 
alysts, sensors or building blocks for molecular electronic 
and spintronic devices [lHa]. They also form an interest- 
ing class of materials where the evolution of the clusters' 
electronic structure and " mctallicity" from molecular to 
colloidal (bulk-like) regimes can be explored, e.g., by 
measuring the optical or electrochemical response with 
varying cluster sizes. 

Recent refinements of the original synthesis^ of MPCs 
by several groups have enabled monodispcrse samples 
of a few particularly stable magic compounds in 1 - 
3 nm range, and a few have by now been determined 
up to molecular precision; these include Au2o(SR)i6 [1], 



Au2s(SR)- 1/0 |4nj, Au 38 (SR) 2 4 HEl, Au 4 o(SR) 24 
[3], Au68(SR)34 15|, Auio2(SR)44 HH, and compounds 
around 144 Au atoms and 60 thiolates [H [H El . The 
AU25, AU38 and AU102 clusters are notable in this series 
since their total atomic structure has been determined 
from X-ray crystallography, opening the door to detailed 
theoretical analysis of the surfacc-covalent gold-sulphur 
bond and the electronic and geometric factors underlying 
the stability of these specific compounds. 

Up to now, these MPCs have been investigated as 
stand-alone entities computationally. However, under- 
standing the interaction of MPCs with their environ- 
ments is fundamentally important in order to assess their 
suitability as building blocks for components in applica- 
tions in molecular electronics or spintronics. Here, we 
have embarked on these studies by studying computa- 
tionally the interaction of the Au2s(SR)i8 cluster and 




FIG. 1. The structure of (a) Au25Lj~ g 1 cluster, (b) AU13 
core, (c) AU2L3 unit and AU25L18 cluster on (d) Au(lll) , (e) 
Ag(lll) and (f) Cu(lll) surface, (g), (h), (i) are the same 
structures as (d), (e), (f), but rotated 90 degrees. The atoms 
are colored as: Au: orange, S: yellow, C: dark grey, H: light 
grey, Ag: silver, Cu: brown. 



its magnetically doped derivative MnAu24(SR)i 8 with a 
prototypical metal surface, for which we select the noble 
metal (111). 

A particularly stable phcnylcthanethiolatc-protected 
gold cluster, Au2 5 (SCH 2 CH 2 Ph) 18 1 was structural! 



characterized from X-ray crystallography in 2008 
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and its electronic and optical properties have been an- 
alyzed by using the density functional theory (DFT) 
fiol fl9l [ioj . The structure of the cluster is shown in 
Fig. QJ, (from here on, we use methylthiolate SCH 3 as a 
model ligand and abbreviate L = SCH 3 ). It comprises 
an approximately icosahedral core of 13 Au atoms (see 
Fig. [Tp) surrounded by six AU2L3 units (see Fig. [Ip) in 
an octahedral arrangement with 12 S-Au contacts that 
passivate the core surface. The cluster Au2sL^ 8 is known 
to be a stable redox species with at least three possi- 
ble charge (and spin) states: q — -1, 0, +1. This has 
led to suggestions that replacement of one or more core 
Au atoms by a magnetic atom such as Mn should lead 
to a "magnetic superatom"; theoretically this has been 
shown to be the case for an isolated MnAu24Li§ 2(J 21 



The robust geometry and a well-understood electronic 
and magnetic structure of these clusters motivated us to 
use Au25Li 8 and MnAu24Li 8 to investigate how metallic 
surfaces will affect the structure and properties of MPCs 
and whether the thiolate ligands can decouple the delo- 
calized superatom orbitals and local spins from the elec- 
tronic structure of the surface. 

We use spin-polarized DFT as implemented in the 
GPAW 



code 
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The DFT calculations arc 
performed with the projector-augmented-wave (PAW) 
method in a real space grid (0.2 A grid spacing) using the 
Perdew-Burke-Ernzerhof (PBE) generalized-gradient ap- 
proximation [H]. We regard Cu(3d 10 4s 1 ), Ag(4d 10 5s 1 ), 
Au(5ci 10 6s 1 ), Mn(3d 5 4s 2 ), S(2s 2 2p 4 ), C(2s 2 2p 2 ) and 
H(ls 1 ) electrons as the valence. The PAW setups for 
all the metals include scalar-relativistic corrections. Pe- 
riodic boundary conditions are applied in all the dimen- 
sions. The Au(lll), Ag(lll) and Cu(lll) substrates are 
modeled with three layers and 15 A vacuum is used to 
separate the AU25 cluster and the image slab in the z 
direction. The two bottom metallic layers are fixed but 
all other atoms are relaxed during the geometry opti- 
mization until forces are below 0.05eV/A. A (5x6) sur- 
face cell is constructed for Au(lll) and Ag(lll) surfaces 
while a (6x8) surface cell is constructed for Cu(lll). A 
Monkhorst-Pack 2x2x1 k-point sampling and the Bader 
method[25| are used to analyze the electronic structure 
and the charge state of the cluster. 

The relaxed structures of Au 2 5Li 8 on Au(lll), 
Ag(lll) and Cu(lll) surfaces are shown in Fig. [TJi-[TJ. 
The geometrical structures of AU25L18 on Au(lll) and 
Ag(lll) surfaces are very similar. The shortest distance 
between the cluster and the surface atoms is 2.66 A for 
the Au(lll) surface and 2.69 A for the Ag(lll) surface 
which is the distance from one S atom to the Au (Ag) sur- 
face. Besides the S atom, there are two Au ligand atoms 
close to Au(lll) and Ag(lll) surfaces at 2.92 A, 2.93 
A from the Au(lll) surface and 2.93 A, 2.98 A from the 
Ag(lll) surface. These three connections stabilize the 
cluster on the surface. Although some SCH3 units are 
rotated by the interaction with the surfaces, the overall 





Cluster 


Core 


Ligands 


Au2 5 L lg 1 


-1.00 


0.31 


-1.31 


Au 25 L 18 /Au(lll) 


0.28 


0.41 


-0.13 


Au25Li 8 /Ag(lll) 


-0.42 


0.34 


-0.76 


Au 25 Li 8 /Cu(lll) 


-0.33 


0.41 


-0.74 



TABLE I. The calculated charge decomposition (by the Bader 
method) of Au2sL7/ g and for the cluster on Au(lll), Ag(lll) 
and Cu(lll). The core consists of the approximatively icosa- 
hedral AU13 and the ligand layer consists of the six AU2L3 
units. 




FIG. 2. The isosurface (h=0.005A~ 3 ) of Sp for Au 25 Li 8 
on (a) Au(lll), (b) Ag(lll) and (c) Cu(lll) surface. The 
blue and red regions show accumulated and depleted electron 
charge, respectively. The colors of atoms are as in Figure 1. 



shape of the metal frame of the cluster remains spherical 
as in the case of Au25Lj~ g 1 in the gas phase. 

The structure of AU25L18 on Cu(lll) is significantly 
different from Au(lll) and Ag(lll). The shortest dis- 
tance between the cluster and the surface is also the 
distance between one S atom and Cu(lll) surface (2.49 
A), but the symmetry of the cluster is severely broken. 
Two ligand Au atoms are displaced towards the surface: 
one stays approximately on the hollow position of the 
Cu(lll) surface and another stays approximately on the 
bridge position. The shorter distance and the distortion 
of the cluster show that AU25L18 and Cu(lll) surface 
have a strong interaction. This is reflected by the calcu- 
lated adsorption energy AE of the cluster on the surface, 
for which we find values of -1.15 eV for Cu(lll), -0.77 eV 
for Ag(lll), and -0.53 eV for Au(lll). It should be noted 
that our values are most likely slight undcrcstimations of 
the true adsorption energies since the PBE functional 
does not inlude van der Waals -type (vdW) interactions. 
Anyway we expect that the above-mentioned systematic 
behavior is not changed by the neglected vdW effects. 

When a MPC is placed on a metallic surface, the in- 
teractions between the MPC and the surface may induce 
the charge transfer effect which might change the elec- 
tronic and magnetic properties of the cluster. In or- 
der to understand this effect, we calculated the Bader 
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charge distribution of the Au 2 5Li 8 cluster on the dif- 
ferent surfaces and compared it to the stable Au25Lj" g 1 
in gas phase. The charge decomposition in the core 
and in the ligand layer of the cluster is given in Ta- 
ble HI 0.28 electrons are transferred from Au2sLi 8 to the 
Au(lll) surface, while 0.42 electrons are transferred from 
the Ag(lll) surface to the cluster and 0.33 electrons are 
transferred from the Cu(lll) surface to the cluster. We 
note that the charge transfer between the cluster and 
the surface shows the same trend as the work function 
(W) of the surface. The order of the work functions 
is W(Ag(Hl))<W(Cu(lll))<W(Au(lll)). If the Bader 
charges of the Au core and ligands are calculated sepa- 
rately, one can find out that the Au core has the simi- 
lar charge state for A^sL^g 1 cluster and Au25L 18 on the 
different surfaces, while the charge states of ligands are 
quite different, which indicates that the charge rearrange- 
ment of AU25L18 mainly happens in the ligand layer when 
the cluster absorbs on the Au(lll), Ag(lll) and Cu(lll) 
surfaces. 

The real-space charge redistribution can be visualized 
by computing the quantity Sp = Pdu+sub - Pdu - Paub, 
where Pdu+sub is the electron charge density for the 
combined system while p c i u and p su f, are the electron 
charge densities of the noninteracting cluster and sub- 
strate, respectively. The isosurface (h=0.005A -3 ) of 6p 
for Au 25 Li 8 on Au(lll), Ag(lll) and Cu(lll) is shown 
in Fig. [2] The blue region is where the electrons are 
accumulated and the red one is where the electrons are 
depleted. It is clear that the charge rearrangement is lo- 
calized mainly in the ligands of the cluster and in the 
first layer of the metal surface. The different adsorp- 
tion mechanism can also be seen in this figure. For 
Au25Lig/Au(lll), the ligands which are close to the sur- 
face lose electrons, while the Au surface atoms which are 
close to ligands gain electrons. The coulomb interaction 
between the positive charge of the cluster and the nega- 
tive charge of the surface makes the cluster stable on the 
surface. For Au25Li 8 /Cu(lll), the electrons accumulate 
between the cluster and the surface which shows forma- 
tion of a chemical-bond-like interaction between the clus- 
ter and the surface. 

It is of great interest to analyze the electronic struc- 
ture of the adsorbed cluster, especially the frontier or- 
bitals near the Fermi level. The Ai^sL^g 1 has been stud- 
ied previously 0, M, 111- The HOMO level is 3-fold 
degenerate whereas the LUMO is 2-fold degenerate with 
a HOMO-LUMO (HL) gap of about 1.3cV (PBE value). 
The character of the electronic states of the cluster can be 
analyzed by projection onto spherical harmonics placed 
at the mass center of the cluster ljl, 2(1 : 




FIG. 3. The angular-momentum-projected local density of 
states (PLDOS) (projection up to I symmetry, i.e I = 6) for 
the AU13 core of (a) Au2sL]~ g , (b) AU25L18 on Au(lll), (c) 
AU25L18 on Ag(lll), and (d) AU25L18 on Cu(lll). The indi- 
vidual weights c iy i have been folded by a Gaussian of 0.1 eV 
width. The Fermi energy is at zero. 



where 



ipi,lm(r) = / dfYi m (r)(j>i(r), 



(2) 
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r 2 dr \ipi,i m (r) 
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Here i is the index of the Kohn-Sham state 4>(r) and Y; m 
is the spherical harmonic function with I as the angular 
quantum number and m as the magnetic quantum num- 
ber. The angular momenta arc considered up to I = 6 
(I-symmetry) and the expansion is made in a sphere of 
radius Ro, where Ro is chosen to be 4.5 A. This analy- 
sis reveals the symmetry of the Kohn-Sham states of the 
system in the Aui3 core. Previously this analysis showed 
for Au25Lj" g 1 that the HOMO has P-symmetry and the 
LUMO has D-symmetry as shown in Fig. [3Jt. This result 
is easy to understand from the superatom model: The 
"free-electron" count for Au 2 5Lj" g 1 is defined by the 26 
Au(6s) electrons (including the extra electron from the 
negative charge) "donated" to the sp-subsystem, out of 
which 18 electrons are withdrawn by the thiolates, result- 
ing in eight delocalized electrons with closed electronic 
shell at lS 2 lP 6 .[l]l[26| 

Here we extend the use of eqs. 1 and 2 to analyze 
the electronic shells of the metal core of the adsorbed 
clusters. The angular- momentum-projected local elec- 
tron density of states (PLDOS) for the AU25L18 on the 
Au(lll), Ag(lll) and Cu(lll) are shown in Fig. H>Hl 
respectively. Neglecting the very small feature just above 
the Fermi level (originating from the interaction between 
the cluster and the surface) , the HOMO of the cluster in 
all cases has a clear P-symmetry and the LUMO has a 
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FIG. 4. (a) The strucure of MnAu 24 Li 8 on Au(lll). Mn 
is green, the other colors as in Figure 1. (b) The angular- 
momentum-projected local density of states (PLDOS) of the 
MnAui2 core (Ro = 4.5A in eq. 1) and Mn atom (i?o = 1-4 A 
in eq. 1) of M11AU24L18 on Au(lll). Top and bottom panels 
show the PLDOS for spin up and down electrons, respectively. 
The individual weights c it i have been folded by a Gaussian of 
0.1 eV width. The Fermi energy is at zero. 



clear D-symmetry. This shows that the " superatomic or- 
bitals" exist also when the cluster is adsorbed on the 
metal surfaces. From the peak positions, the appar- 
ent "HL gaps" for the clusters can be determined. For 
the Au(lll) and Ag(lll) surfaces, the HL gap is about 
1.0 eV, slightly reduced from the gas phase value of the 
Au25LJ~g 1 cluster. The HL gap further diminishes when 
the cluster is on the Cu(lll) surface (0.7eV), mainly 
because the LUMO states broaden into a wider band, 
caused by the stronger interaction between the cluster 
and Cu surface. 

If a cluster retains a magnetic moment on the sur- 
face, it may have possible applications for spin-dependent 
transport in molecular spintronics. The neutral AU25L18 
has a magnetic moment 1 fis in the gas phase, but we 
found that the moment is quenched on all the (111) 
surfaces. The Mn-doped MnAu24Lig has a high mo- 
ment of 5 up in gas phase due to the localized Mn(3d 5 ) 
electrons. 20|, 21 1 We found that the energetically opti- 
mal site of the Mn atom is in the surface of the Aui3 
core, in agreement with another recent DFT study (2 71]. 
The energy differences to the next structure isomers are 
0.30 eV and 0.36 eV (corresponding to Mn in the ligand 
layer and at the center of the core, respectively). Mo- 
tivated by this result we replaced one of the Aui 3 core 
surface atoms, facing Au(lll), by the Mn atom. We find 
that the structure of the doped MnAu24Lig/Au(lll) and 
Au25Li8/Au(lll) are very similar (see fig. H^). However 
in contrast to the Au25Lig/Au(lll) case, we observe that 
the magnetic moment of MnAu24Lis/Au(lll) remains 
very high, 4.9 fiB, with 4.6 [is localized on the Mn atom. 
We show the PLDOS of the MnAu i2 core and the Mn 
atom of MnAu2 4 Li 8 /Au(lll) in Fig. 0Jd. The 3d or- 
bitals of the Mn atom remain fully spin-polarized which 
gives rise to the high magnetic moment of the cluster. 



Simultaneously, the frontier orbitals of the cluster retain 
the P symmetry for HOMO and D symmetry for LUMO. 

In conclusion, using mcthylthiolatc protected AU25L18 
and a doped MnAu24Lis as prototype clusters, we stud- 
ied the interaction between MPCs and noble metal (111) 
surfaces. The geometric structure of the cluster and the 
adsorption mechanism depend on the interaction between 
the ligands and the surface, but the structure and the 
electronic properties of the core remain. The charge rear- 
rangement is localized in the ligand layer which protects 
the overall electronic shell structure of both AU25L18 and 
MnAu24Lig and localized spins of Mn 3d electrons in the 
core of MnAu24Li 8 by decoupling them from the surface. 
The HOMO-LUMO gap of the adsorbed cluster depends 
slightly on the metal, being the smallest on the Cu(lll) 
which also has the strongest interaction with the cluster. 
However a clear gap of about 0.7 eV to 1 eV is observed 
in all cases. Our results implicate that monodispersed 
thiolate-protected gold clusters retain their "superatom 
identity" also when adsorbed on prototypical metal sur- 
faces. These novel building blocks can be experimentally 
synthesized in pure forms, can be doped by magnetic 
atoms, and their ligand layer can be functionalized by 
varying the nature (e.g., aryl- or alkylthiolates) and size 
of the ligands which affects the electronic coupling and 
tunneling barriers. Utilizing this diversity may open new 
ways to engineer building blocks for molecular electronic 
and spintronic devices. 
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